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We present a measurement of the longitudinal spin asymmetry A\\ in photoproduction of pairs of 
hadrons with high transverse momentum pr- Data were accumulated by the HERMES experiment 
using a 27.5 GeV polarized positron beam and a polarized hydrogen target internal to the HERA 
storage ring. For h + h~ pairs with p^ 1 > 1.5 GeV/c and p 1 ^? > 1.0 GeV/c, the measured asymmetry 
is An = —0.28 ± 0.12 (stat.) ± 0.02 (syst.). This negative value is in contrast to the positive 
asymmetries typically measured in deep inelastic scattering from protons, and is interpreted to arise 
from a positive gluon polarization. 



From polarized deep inelastic lepton-nucleon scatter- 
ing (DIS) experiments, it has been inferred that the quark 
spins account for only a fraction of the nucleon spin. One 
possible explanation is a significant gluon polarization in 
the nucleon. In principle, the polarized gluon distribution 
AG(xg) {xq is the fraction of the nucleon momentum 
carried by the struck gluon) can be probed by a mea- 
surement of the scaling violation of the polarized struc- 
ture functions. However, the presently available data on 
polarized inclusive deep inelastic scattering only poorly 
constrain AG(xg), although there is some indication for 
the integral to be positive . On the other hand, two 
theoretical calculations in the bag model obtain different 
predictions for the sign of the integral of AG(x G ) §§• 
Several recent experimental proposals have concentrated 
on ways to measure AG(xg) directly |§-|J. 

One way to measure AG(xg) directly is via the photon 
gluon fusion process. Two useful experimental signatures 
of this process are charm production and production of 
jets with high transverse momentum px- In the former 
case, the large mass of the charm quark suppresses its 
production in the fragmentation process. A similar argu- 
ment applies to the production of jets: the transverse mo- 
mentum produced in the fragmentation process is small 
and two back-to-back jets with sufficiently high px thus 
reflect the high pt of the quark and anti-quark produced 
in the photon gluon fusion process. Both charm produc- 
tion and high-pT jet production in DIS have resulted in 
direct measurements of the unpolarized gluon structure 
function G(x G ) §-f|. 

At lower energy fixed target experiments, high-p^ 
hadrons must serve in place of jets [^2| . Several phe- 
nomenological studies of the potential of high-p^ meson 
photoproduction as a probe of AG(xg) have been per- 
formed HQ. 

In this Letter we present the first measurement of 
a spin asymmetry in photoproduction of pairs of high- 
Px hadrons. The data were collected in 1996 and 1997 
by the HERMES experiment at the HERA storage ring 
of the DESY laboratory. Polarized positrons of energy 
27.5 GeV were scattered off a polarized internal hydro- 
gen gas target. The beam polarization was continuously 
measured by Compton back scattering and had an aver- 



age value of 0.55 ±0.02 pj|jl6|. The average target polar- 
ization was 0.86 ± 0.04 ]15| , |17| . In both cases the quoted 
uncertainty is predominantly systematic. The HERMES 
detector jig] is a forward spectrometer that identifies 
charged particles in the scattering angle range of 0.04 
< 8 < 0.22 rad. Particle identification (PID) is accom- 
plished using an electromagnetic calorimeter, a scintilla- 
tor hodoscope preceded by two radiation lengths of lead, 
a transition radiation detector, and a C4F10/N2 (70:30) 
gas threshold Cerenkov counter. A likelihood method, 
based on the empirical responses of each of the four 
PID detectors, is used to discriminate between positrons 
and hadrons. The luminosity is measured in a pair of 
NaBi(W04)2 electromagnetic calorimeters that detect 
Bhabha-scattering from target electrons. 

The longitudinal cross section asymmetry A\\ was de- 
termined using the formula: 

Au = -f rr- . (1) 

Here TV^(jV^) is the number of oppositely charged 
hadron pairs observed for target spin parallel (anti- 
parallel) to the beam spin orientation. The luminosities 
for each target spin state are L^^^ and Lp^^\ the lat- 
ter being weighted by the product of the beam and target 
polarization values for each spin state. 

Events were selected that contained at least one pos- 
itively charged hadron h + and at least one negatively 
charged hadron hr . The observation of the scattered 
positron was not required, in order to include the very low 
Q 2 region which dominates the cross section. Here Q 2 is 
the negative square of the 4-momentum of the virtual 
photon. The highest momentum hadrons of each charge 
were required to have a momentum above 4.5 GeV/c and 
a transverse momentum px above 0.5 GeV/c. Here px 
is defined as the momentum transverse to the positron 
beam direction and is approximately equal to the mo- 
mentum transverse to the photon direction when Q 2 « 0. 
To suppress contributions from vector meson resonances 
from the data sample, a minimum value of the invariant 
mass of the two hadrons (assuming both hadrons to be 
pions) M(2-7r) > 1.0 GeV/c 2 was imposed. Additionally, 
both hadrons were required to originate from the target 
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region and to have a common vertex. A detailed account 
of the analysis may be found in |19] . 
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Fig. [l] presents the measured A\\ for the highest values 
of transverse momenta accessible at HERMES; in the 
top (bottom) panel the positive (negative) hadron was 
required to have a pr greater than 1.5 GeV/c and the 
asymmetry A\\ is then plotted as a function of the px 
of the hadron of opposite charge. The data suggest a 
more negative asymmetry when the transverse momen- 
tum of the negative hadron is higher than that of the 
positive hadron. Ignoring this charge asymmetry and 
averaging over the five bins satisfying the requirement 
p't > 1.5 GeV/c and p^? > 1.0 GeV/c, a negative asym- 
metry A| = -0.28±0.12 (stat.)±0.02 (syst.) is observed. 
(The symbol hi signifies the hadron with the higher pt-) 
When the requirement p^ > 1.5 GeV/c is not enforced, 
the asymmetry is consistent with zero. The observed 
negative asymmetry is in contrast to the positive asym- 
metries typically measured in deep inelastic scattering 
from protons. 

A possible background to the observed asymmetry 
arises from coincident detection of a negative hadron 
and the scattered positron, the latter being misidentificd 
as a positive hadron. From studies of other processes, 
the probability for positron/hadron misidentification has 
been determined to be less than 0.2%. By comparing 
yields of h + h~ pairs to those of e + h~ pairs detected in 
the final state, the background arising from this misiden- 
tification has been estimated to be less than 0.1%, for 
the kinematics selected by this analysis. Other sources 
of background include high-py particles from charm de- 
cays. Contributions from both open charm and J /if) de- 
cays have been found to be negligible using the AROMA 
p0{ Monte Carlo generator. 

The systematic uncertainty arising from the measure- 
ment of the beam and target polarizations is about 6% of 



Am , much smaller than the statistical error and indepen- 
dent of pt- Resolution effects and alignment uncertain- 
ties were found to be negligible. Electroweak radiative 
corrections are expected to be very small compared to 
the statistical uncertainty. 

The measured asymmetry was interpreted assuming 
that several different processes could contribute to the 
two-hadron cross section: lowest order deep inelastic 
scattering (containing no hard QCD vertex), interaction 
via the hadronic structure of the photon - described 
by the vector meson dominance model (VMD) and by 
non-resonant hadronic "anomalous" photon structure, 
and the two first order QCD processes (termed "di- 
rect") which describe the interaction of a pointlike pho- 
ton. These are photon gluon fusion (PGF) and the QCD 
Compton effect (QCDC). 

The contribution from lowest order DIS is suppressed 
by the requirement of high pr, and was confirmed to be 
negligible by a simulation based on the LEPTO Monte 
Carlo generator |2lj . Contributions from VMD were as- 
sumed to have a negligible spin asymmetry, and were 
thus treated as a dilution of the other asymmetries. Fi- 
nally, we neglect possible contributions from anomalous 
photon structure, where the photon fluctuates into a non- 
resonant qq pair which interacts via hard processes with 
the partons inside the nucleoli. This is supported by a 
model p3] that explains the excess of forward hadrons 
with high pt observed in jp reactions at 70 — 90 GeV, 
relative to those from irp and Kp scattering |23| . At this 
energy, the model prediction at high pt is dominated 
by direct processes involving hard coupling of the pho- 
ton to the partons in the proton. At the lower energy 
of the present experiment, a negligible contribution from 
anomalous photon structure is predicted by the model. 

Under the assumptions described above, only two of 
the five possible spin asymmetries A\ contribute signifi- 
cantly to the measured asymmetry: 



A\\ fa (Apgf/pgf + Aqcdc/qcdc)-D 



(2) 



where /, is the unpolarized fraction of events from sub- 
process i (/pgf + /qcdc + /vmd = 1), and D is the 
virtual photon depolarization factor. In the small region 
of phase space selected by the present analysis, the A^s 
may be approximated by the products of the hard sub- 
process asymmetries and the quark and gluon polariza- 
tions. The subprocess asymmetries opgf = a(7<7 - > (?<?) 
and q.qcdc — o,{lQ ~ * 19) arc directly calculable in lead- 
ing order (LO) QCD [^3|. For real photons and massless 
quarks, Spgf = — lj while (aQCDc) is about +0.5 (aver- 
aged over the kinematics selected by this analysis) and 
is independent of the quark flavor. The effective quark 
polarization Aq/q is computed as a suitably weighted 
combination of Au/u and Ad/d, which are known from 
inclusive and semi-inclusive polarized DIS measurements 
p4| , p5| . The measured asymmetry can therefore be ex- 
pressed as follows: 
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A 



|| ~ ^OPGF-^-/pGF + aQCDC-^/QCDC 



where the kinematic dependences have been suppressed 
for brevity. This equation can be solved for AG/G after 
appropriate averaging over the selected kinematics. 

The PYTHIA Monte Carlo generator [|(| was used to 
provide a model for the data. An important parameter 
in the simulation of the direct processes, the minimum 
transverse momentum of the outgoing partons (pjp m ), 
was chosen following Ref. p2j to be 0.5 GcV/c. The 
kinematic region used in the interpretation of the mea- 
surement (p^ > 1.5 GcV/c and p%? > 0.8 GcV/c) was 
chosen so that the final results depend only weakly on the 
choice of p™ m . The Lund fragmentation parameters used 
in the simulation have been adjusted to fit the HERMES 
semi- inclusive hadron multiplicity data fl27j, 




FIG. 2. Comparison of data (circles) and Monte Carlo 
simulation (full histogram) for dN/dp? 2 for p^, 1 > 1.5 GeV/c. 
The dashed, dashed-dotted and dotted lines represent the con- 
tributions from the PGF, VMD and QCDC processes, respec- 
tively; the solid line represents their sum. 

The normalized yield for the production of two high- 
Pt hadrons is compared to the Monte Carlo simulation 
in Fig. ||. Here, the Weizsacker- Williams approximation 
has been used to relate the photoproduction cross section 
simulated by PYTHIA to the measured clcctroproduc- 
tion cross section. Also shown in Fig. || are the contribu- 
tions from the three subprocesses included in the simu- 
lation. The simulated yield has a pr-dependence similar 
to that of the data, but is significantly smaller in mag- 
nitude. Good agreement is found for the distributions 
in other kinematic variables, such as the azimuthal angle 
between the two hadrons and Apt = \p^ \ — |py |. As 
the simulation of the direct QCD processes is restricted 
to leading order, the observed difference in normalization 
might be due to contributions from higher-order QCD 
processes and/or contributions from hard interactions of 
the hadronic structure of the photon. We also note that 
the agreement becomes much better if the default Lund 



fragmentation parameters are used. However, the final 
result for AG/G is found to depend only weakly on the 
choice of fragmentation parameters. 

In the same region of phase space where a negative 
asymmetry is observed (p^ 1 > 1.5 GcV/c and p^j? > 
1.0 GeV/c), the simulated cross section is dominated by 
photon gluon fusion. The consequent sensitivity of the 
measured asymmetry to the polarized gluon distribution 
is demonstrated in Fig. |^, where A| at high transverse 
momenta (i.e. the average of the two panels of Fig. |l|) 
is compared with Monte Carlo predictions for different 
distributions of AG/G. 
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FIG. 3. An for high-pT hadron production measured 
at HERMES compared with Monte Carlo predictions for 
AG/G = ±1 (lower/upper solid curves), AG/G = (mid- 
dle solid curve), and the phenomenological LO QCD fits of 
Ref. psj (dashed, dotted, and dot-dashed curves). 



The Monte Carlo simulation was used to determine 
the quantities necessary to relate the data to (AG/G). 
where the angle brackets indicate averaging over the kine- 
matics of the measurement. These quantities have been 
determined from PYTHIA to be (Daqcnc^) = 0.15, 

(D) = 0.93, (x G ) = 0.17, (Q 2 ) = 0.06 (GeV/c) 2 , and 
ip%) = 2.1 (GcV/c) 2 . The distribution AG(x G ) is 
probed principally in the range 0.06 < xq < 0.28. Note 
that the hard scale of this process is not given by Q 2 , 
but rather by p\, the square of the transverse momen- 
tum carried by each of the outgoing quarks. 

For the four values of A\\ at Pj? > 0.8 GeV/c presented 
m Fig. | (AG/G) was extracted according to equation 
(^). Since these four measurements probed essentially 
the same range of xq, the results for (AG/G) were av- 
eraged. Using the assumptions and model parameters 
described above, a value for (AG/G) was determined in 
LO QCD to be 0.41 ± 0.18 (stat.) ± 0.03 (syst.), where 
the systematic uncertainty represents the experimental 
contribution only. 

The extracted value of (AG/G) is compared in Fig. ^ 
with several phenomenological LO QCD fits of a subset 
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of the world's data on gi(x, Q 2 ) p8| , p9[ . The horizon- 
tal error bar represents the standard deviation of the xq 
distribution for the cited kinematical constraints on the 
produced hadrons, as given by the Monte Carlo. 
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FIG. 4. The extracted value of AG/G compared with 
phenomenological QCD fits to a subset of the world's data on 
gi' n (x, Q 2 ). The curves are from Refs. evaluated at a 

scale of 2 (GeV/c) 2 . The indicated error on AG/G represents 
statistical and experimental systematic uncertainties only; no 
theoretical uncertainty is included. 



In summary, a positive value for the gluon polarization 
has been extracted from a measurement of the spin asym- 
metry in the photoproduction of pairs of hadrons at high 
Pt- This interpretation of the observed negative asym- 
metry is based on a model which takes into account lead- 
ing order QCD processes and VMD contributions to the 
cross section. At the kinematics of this measurement, no 
spin-dependent analyses of higher order QCD processes 
or contributions from anomalous photon structure are 
presently available; these processes have therefore been 
neglected in the model presented here. If such processes 
would be important but have no significant spin asymme- 
try, the extracted value of (AG/G) would increase, but 
still differ from zero by 2.3cr. To alter the principal con- 
clusion of this analysis, i.e., that (AG/G) at (xq) = 0.17 
is positive, a significant contribution from a neglected 
process with a large negative spin asymmetry would be 
needed. 
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funding agencies for financial support. 
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